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ABSTRACT 


The successful application of holographic interferometry to the 
study of density fields around opaque bodies in wind tunnel experiments 
has been reported in the literature. The present report extends this 
technique to the study of the three-dimensional asymmetric flow fields 
encountered near the wing-fuselage junction of an aerodynamic model 
in the transonic flow regime. Finite fringe interferometry has been 
used to obtain fringe information about a partially transparent wing- 
body structure. A FORTRAN computer program was utilized to invert 
the fringe information and produce a plot of the density field around 
the model. The resulting asymmetric density field and shock wave 
structure are shown to be an accurate representation of the phenomena 


encountered in aerodynamic corner flow. 


Nima? 
‘gorse afeuberate14 tives 
QPQES simatiied you all 


. “1? af 29 @edpaes oat 


Hes Lanne Wee ee 


tet} ‘éa Jteqe- aoline aaa ind ae 
pf ': #h yWesptrnsgiea aot tna PSA ails Yo 


weer ee 1:70 tana 
. Ce etnies 


| st a 


tees theltym «tie se ets opi 


eis oa ig na OE A or 
Al ate 61) meal 9 Naan eae ed 


9 


ie Wy) eteemaae. Pie Larter en palttonsdl 


16 oO; dome eeneeinaalan wid: 


. ~~ “ . an 
* WOR sherta gctaleghat eee 


=? 


EL. 


Hit, 


Vs 


Vig 


Vas 


TABLE OF CONTENTS 


LIN CIENE ON DNOM SARIN O INS OE IS cS a gg eee te a a a a 
pe Parolnyeb Nig teal PARA US) eis 2 54 Sse as Sek 
Peed e Ly UD eee NING = <i  eee  Do o ee S 
pote IOrmOGk A PHICsARRANGEMENT -osLe..o--5-.2- 
Carlee WIND Me NNELUMODE Lice -s522ssceueck ao sel 
ANAT Y LIGA EVALUATION*‘OFP THE DENSITY FIELD-2-- 
Awe BASIC EQUATION OFINTERPEROMETRY ------ 
ee ree eee Gr Eeelet NIV RoIONS S22 ot See oa bssese So Se 
Ca heen iCAat PROCEDURE et 22c: J2ces. Ss Soe 
Pere VUE aden NOG DUI Ba aS eee Sec oe 
aD OLA OK Ye Ll WC LNIOUMS = =o atm ete ec ees 
Bee HOmOghweAre MiG. TmCHNIOQUES |. -22---oscccnsese een 
Cr wie OCLC nM OIN aacre sae aoa wea aoe oes nok es Sie 
Peep RIMEN TAL RESULTS AND DISCUSSION -2.=-+-.5.-< 


INCAS ON Se 


APPENDIX A: Reduction Of An Interferogram To Obtain 


Brin cgeqedintWata= aaj Se saae alee aia ofS ok 


ALP ENDIX B: Application Of Computer Program 'HOLOFER'- 


SLA SM IC CQUISE 18U1E5 13 D520 SIN Cs se eo 


PiiVMee ee Miomii OM ON MLS: j. fs saonacea seen esa ee sees tes 


TET ONS Sal TEND U2 o/c 


i L 


Bie 


_ 


7 4 : ¢@ 
ile en ay = ib ai arte Pf 
. g 7 
oi ’ - 
yp Pi 7 ~ SS. : > 
7 7 : 
. 2a he TPA Ee bai cb j . 
‘ Pree "ww tle roay iA iD LOuE tre 
: mae ' oe ‘- 7 | 
‘ ; o) eal ke We a 
a - oe e\6 @ whe Cie) ee b 
* 1 - =) 
i mn | 
, Ber: Ct Git a worry 
: : 
7 ’ 7 
, req 24) *] ra 
m - 
a er ' ‘ 
~ « - . 5] ‘ 
ine do a e 
, ibaa 
: = oN ity ir ae 


4 _ 
: a . a - 
Rte a » 


ii lide — . Srcay 7: = _ af a 
a 4 4 Cae, I OAD | 
Py ; Te i r ey 
OTOn 
‘ 
4 * ; 
= a 


ACKNOWLEDGEMENTS 


This writer wishes to gratefully acknowledge Dr. D. J. Collins 
for his most valuable guidance and assistance during the course of this 
investigation; the Aerodynamics Division of the Naval Ship Research 
and Development Center, Carderock, Maryland, under Dr. R. Furey, 
particularly Baeihstand=©-., Gilmore, for their tireless:effort during 
the wind tunnel testing phase; the technical staff of the Department of 
Aeronautics under R. Besel and T. Dunton, particularly N. Leckenby 
and G. Middleton, for their able assistance in the fabrication and 
assembly of the models and equipment used in this study; and my wife 


and family for their noble patience and encouragement. 


ee ie W va a | ui 
"i i Vy, 
Hh 0 OE RE ED Ques iy rd Bia ian bb Pa aidgif 


See 
pis SS 1 ew lee iwi won tt ivbt win mid ait a 
i a0 stirtsg - bedi evnaitite Aeron - obi oe 
o nite szaletew te ‘aa +h inpoetll iD bein ait | 
yi 0) tt Le Fie Lei al taney vue 


vir noe rey va vei ie imi. mer oft 


levied ih soe Aaledime hip Fer vat 
f. 
rbiye in) evo a itn bags a" Se te 
. ry 


AP pi! ti iS nie aia rome Ling ad ‘net % ays 


I. INTRODUCTION 


The field of flow measurement has been revolutionized in recent 
years with the perfection of holography and holographic interferometry 
techniques. High power Q-switched and dye-switched lasers and 
sophisticated double-pulsing trigger mechanisms provide exposure 
times on the order of twenty nanoseconds, thereby ''freezing'' the flow 
during the eae production process. The precision optical quality 
components and measurement techniques of Mach-Zehnder inter- 
ferometry have given way to the much less restrictive requirements of 
holographic interferometry which provide high quality interferograms 
in three dimensions. 

Techniques for the application of holography to interferometry 
have been reported by Heflinger, et al. [1], and by Brooks [2]. In the 
determination of the density field around a free jet in the supersonic 
regime, Matulka [3,4] expressed the fringe data ina series of ortho- 
gonal polynomials and transformed them to polynomials representing 
density using an inversion technique reported in [5,6]. The method 
was extended by Jagota [7, 8] to the determination of the three- 
dimensional density field around a ten-degree half angle cone ina 
supersonic wind tunnel. The ability to produce readable holograms 
in wind tunnel studies using transparent phase objects was verified 


by Heyer [9]. Inthe present report the aforementioned techniques 


Lf Ite AT dary ; 
> | i is oe ae : it 7 
att i ay 
a a 
" ih nei 
a 5 ¥ i 
ye et oe) fs sae 
Yao a) Seafeoad ul ihandl 
| ns 7 a ‘ tees 
wOCeTTe Shia peli tae yiqz Wit -2js yu) sag) 
a 7 — im i- 
ruin iwe.@eh Dae Susagle-O iow ay dyth 
2 ; : : ¢ 
¥ jig & ae i y | ry! : . Ben al 
, +i oe sits Atte “gepemllball lard mi vay dels 
Er “ony ggoyt! aa ey eat » eh pn ee WP tes * . 
go 9? ty i a a ir 
si a % a ae bob hoa 
' i omeelti th So% Baheberc ood 
bed eed (eo eM 1 ao pateet Dede ooae conte 
* 


‘| ‘ ; 
ee rr a 4854) vom coVES at 
. Us 


4 rag “VO V8 shila } edocs +@ zi 


a -) 
aed 
TST Te Fre Te Ee se es 
~*~ ’ 5 
i | Psi ‘ P vi ae | ifs if 
1 Fs! { ab; 
ifs CA KGS ¥ 
~ 
th hy { we ty = 


tl} A itil , 


wed yr Oe iti ie | 

1% a iatony 
: a a" i 
' 


’ , 
eee iho | triste ce - HH 
a ad rn SS 


fi . ; | 


Bo i, a ae 


have been combined to study the three-dimensional density field near 
the wing-fuselage junction of an aerodynamic model in transonic flow. 
The experiment was conducted at the Naval Ship Research and Develop- 
ment Center, Carderock, Maryland in an eighteen inch transonic blow- 
down wind tunnel at a Mach number of 0.937, using a semi-transparent 
model of original design. 

Since reasonably small variations in density were anticipated, a 
finite fringe technique wae eee in obtaining the interferograms. The 
horizontal finite fringe field was produced by a vertical translation of 
a diffusing glass in the scene beam a distance of 0.003 inches between 
the two exposures of the holographic plate. Fringe data obtained from 
the interferograms were reduced to density information using a mod- 
ified form of the inversion computer program used in[7]. A self- 
testing procedure incorporated in the program verified the resulting 
density data as an accurate representation of the actual flow around 


the model. 
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ii. EXPERIMENTAL APPARATUS 


A. THE WIND TUNNEL 

The investigation was conducted in the Naval Ship Research and 
Development Center blowdown supersonic wind tunnel. Transonic 
flow conditions were produced through incorporation of slotted upper 
and lower tunnel walls. The tunnel is a nominal eighteen inch blow- 
eto: yacuum facility with a test section fourteen inches by eighteen 
inches in cross-section and twenty-nine inches in length with the slotted 
surfaces installed. Optical quality windows twenty-two inches in 
diameter in the side walls provided complete viewing of the flow in 
the a section as the model was rolled through 180 degrees for holo- 
gram production. A functional schematic of the wind tunnel is shown 
an Figure 1. 
B. THE MACH NUMBER AND PRESSURE MEASUREMENT 

PROCEDURE 

The Mach number at the test section is determined as a function 
of total and static pressure measurements and is maintained by 
carefully controlled butterfly valve settings. Total pressure is 
determined by recording atmospheric pressure prior to tunnel 
operation and accounting for the pressure loss between the plenum 
and the test section during operation. Static pressure is measured 


directly at a central wall port in the test section. Data recordings 
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were made ona Beckman Instruments Company 210 Digital Recorder, 
shown in Figure 2, and were read out on line via a Franklin Strip Tape 


Printer. 


Go LHE HOLOGRAPHIC ARRANGEMENT 

The holographic arrangement is illustrated in Figure 3 and shown 
in photographs included as Figures 4, 5, 6, and 7. Two large wooden 
tables were constructed and linked together with two-by-four beams 
under the tunnel to form the experimental platform. Thick rubber 
pads were attached to the table legs to dampen possible floor vibrations. 
The bulk of the platform provided sufficient stability and vibration 
damping for the experiment. The monochromatic light source used 
was a KORAD K-1 pulsed ruby laser operating at a wavelength of 
6943 Angstroms, together with a Pockels cell Q-switching device. 
The resultant effective exposure time was approximately twenty nano- 
seconds, eliminating the problems due to possible model vibration 
during hologram exposure. To maintain the laser head and output 
etalon at a constant temperature of 27.0 degrees Centigrade, a 
LAUDA constant temperature circulator Model K2R was used. 

The reference beam was directed under the wind tunnel by four 
front surface mirrors, and the beam size was manipulated by means 
of lenses (Figure 3). The scene beam was routed through the test 
section to intersect the reference beam on the holographic plate at 


an angle of approximately fifty degrees. A diffuse glass, mounted on 
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a precision X-Y translation table in the scene beam, was used to 
produce light field holograms. Alignment of the Q-switched laser and 
system optics was accomplished using a continuous wave, low-power 
helium-neon laser. Reference grids were mounted accurately on the 
outer surfaces of the tunnel windows using a survey's transit. Details 
of the model mounting and reference grids are shownin Figure 7. To 
enable hologram production during daylight hours, the entire tunnel 


room was blacked out using drop curtains and light shields. 


D. THE WIND TUNNEL MODEL 

The aerodynamic model used is shown in Figures 8, 9, and 10. 
The metal portions of the model were stainless steel. The greater 
part of the modified double wedge platform wing was constructed of 
optical lucite, as was the portion of the fuselage at the wing root. 
Detailed model dimensions are shown in Figure 1l. The choice of 
aerodynamic design provided good flow characteristics and a relatively 
stable lambda-type shock wave on the wing; the largely transparent 
construction facilitated hologram production through 180 degrees of 
view. 

The model was rotated about its sting mount in the wind tunnel 
from the zero degree position, wings level, to the 180 degree position, 
wings level inverted. Alignment for the desired rotation angle was 
accomplished by manually aligning prescribed lines on the sting mount 


collar with a scribed mark on the sting support. 
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ii ANALYTICAL EVALUATION OF THE DENSITY FIELD 


A. THE BASIC EQUATION OF INTERFEROMETRY 
Interferograms are created when two originally coherent light 
beams are superimposed and projected ona viewing screen. The two 
rays will reinforce or annul each other, depending on their relative 
phase difference at the screen. This phase difference is directly a 
function of the optical pathlengths traversed by the two waves. 
Consider a coherent beam which is split and then recombined ona 
viewing screen. A difference in optical pathlengths of the two com- 
ponent beams may be achieved by causing the beams to traverse 
through different media prior to recombination, with their physical 
pathlengths maintained equal. Each component beam will travel ata 


speed Cy/n where c, is the speed of light ina vacuum and n is the 


0 
index of refraction of the medium traversed. The difference in optical 
pathlength is then given by 

L=L(n, - nj) =c)At (1) 
where At is the time difference of travel in the two media. If the 
optical pathlength is changed by an amount NA , where 7“ is the 
wavelength of the light source and N is an integer, then the order of 
interference changes by an amount N. In other words, a shift of N 


fringes occurs in the interference pattern. The fringe shift may be 


expressed as 
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= L/A (2) 
where g = fringe shift 
“A= light source wavelength 
L= change in optical pathlength 
Substituting equation (1) into equation (2) yields 


g= L (n, - ny) | (3) 


A 


The index of refraction for a given medium is a function of 
density. In the case of gases, since the speed of light is very nearly 
the same as ina vacuum, the index of refraction is well represented 


by the first two terms of a Taylor series expansion [10]: 


n= 14-5 j (4) 


where B dimensionless constant related to the Gladstone-Dale 


constant by K = P/9. 
e = reference density at 0° C, 760 mm. Hg. 
The value of p for air at 7\ = 5893 Angstroms (deep red light) is 
0.000292; variation with wavelength is very small. 
For a fixed difference in the index of refraction between the two 


component beams the fringe shift relation becomes: 


zg = B ae (5) 
S 


For variable density in the test section, the net change in optical 


pathlength is the integrated effect along the beam path, or 


8 L 
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where: Q= 2) Wes (7) 
ance 
fez) Oy (eayazel = 1 (8) 
co 
Z. = plane of constant z 


ds = incremental distance along beam path 
Equation (6) is the basic integral equation for the unknown density. 
With known fringe shift values from an interferogram, the equation is 


inverted to obtain the density along a beam path. 


B. THE INTEGRAL INVERSION 

The integral inversion procedure utilized in this investigation 
was first reported by C. D. Maldonado, et al[5,6]. It was used 
subsequently by R. D. Matulka [3] and R. C. Jagota [7] to determine 
the density variation in an asymmetric free jet and about a cone at 
angle of attack in supersonic flow, respectively. The procedure 
involves the representation of the function f(x,y, Zo) of Equation (6) 
in a complete set of orthogonal functions, with the expansion coeffi- 
oAgHie evaluated by use of the orthogonality condition between the 
functions f and g of Equation (6). The set of functions used is ortho- 
gonal over the entire (x,y) plane for every Zs and remains an ortho- 
gonal set under a rotation of the coordinate system. The coordinate 
system used for the inversion is shown in Figure 12. It consists of 


(a) a set of fixed coordinates x,y and (b) a set of moving coordinates 
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x', y' in which the fringe number function is defined and which rotates 
with respect to x,y as the viewing angle through the test section is 
varied. 
In operator form, Equation (6) can be represented as 
' = 

BA Grey 'yz isl tty, 2.) (9) 

and f is evaluated by inversely transforming the equation to obtain 
= 

Eecnys zy = inet 'e 1,45 Ze) (10) 

This inversion is achieved by utilizing a pair of orthogonal polynomials 
an9-a 


ISAS ee (x, xy) and H,.42;( y') which are related by the trans - 


form relationship 
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where Eo o.y') are Hermite polynomials of order mt+2k. The 


unknown function f(x, y, Zs) is expanded ina set of functions Ween 
as 
fuyny= > F eNO 
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sal i = = 
where E = i for m = 0, ee: Petoreme=. is 2. 3) 6. .,nand 
oe are the unknown coefficients of the expansion, Ct is an 
arbitrary scale factor which may be considered the reciprocal ofa 


non-dimensionalizing coefficient. 
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Utilizing the transform relation of Equation (11), Equation (6) 


becomes 
ye + 1 
4(5,Y, 2) = Vhs a Em| k! (m+k)| gem ste /2, 
mzo K=O 
uy - = ~ofty'2 (13) 
° cms (x) aus oe (x) ans | eee (x) e 
Equation (13) is subject to the orthogonality eoudution 
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where df is the Kronecker delta function. The solution of Equation 


(14) applied to Equation (13) yields the series coefficients 


“17 @ 
tm +k) | ime, 
=e (tlle | (a(Gyi2H, eve ydg as 
“Ny 


With the substitution of the coefficients of Equation (15), Equation (7) 
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The functions Uli are defined as 
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Equations (19), (20), and (21) are the basic equations used to obtain 
the density distribution from the experimentally determined fringe 


variations in a completely asymmetric flow field. 


€. THE NUMERICAL PROCEDURE 

Because the function g( 5 ewes z is an experimentally deter- 
mined quantity the unknown coefficients Bro | ot ) and Ie Of) 
in the series representation of f(x,y, z.) in Equation (19) cannot be 
calculated analytically. It is therefore necessary to evaluate the 
' double integrals of Equations (20) and (21) numerically. This is 
accomplished by noting in Figure 12 and Equation (8) that there is an 


area outside which the density is invariant, namely outside the test 
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section where the known density is Pp or, Since the function f(x, \y;, Ze) 
= 0 outside this circular domain, the limits of integration of + o and 

- oo in Equations (20) and (21) can be replaced by finite values. The 
fringe distribution is then approximated by small increments over the 
test domain, resulting in the representation of the B and D coefficients 


as double series: 
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Using the derivative formula for Hermite polynomials, Equations (22) 


and (23) can be manipulated to yield workable series expressions: 
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Since it is impossible to sum over an infinite number of terms, 


(24) 


(25) 


Equation (19) is necessarily expressed as the sum ofa finite series: 
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It has been demonstrated that judicious selection of the parameters 
Az , Mx, K, M, and & yields density distributions with very 


good accuracy [3, 6]. 
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IV. EXPERIMENTAL PROCEDURE 


A. LABORATORY TECHNIQUES 

In order to visualize the general flow patterns and localize shock 
Or expansion waves about the modela series of standard Schlieren 
photographs were taken at varying flow Mach numbers. Pictures were 
produced for roll angles of ‘Chee 45°, and 90° at Mach numbers from 
0.925 to 1.10. A representative series of Schlieren photographs is 
shown in Figures 13, 14, and 15. Analysis of the Schlieren photo- 
graphs dictated a flow Mach number of 0.937 for the experimental 
study; this Mach number yielded uniform upstream flow conditions and 
located the’ lambda-type shock wave ideally near the center of the lucite 
section of the model wing. 

The coherence length of the pulsed ruby laser was approximately 
ten centimeters for the output power utilized. ‘his reduced the nor-— 
mally critical requirement for pathlength equality in the scene and 
reference beams that must be fulfilled in ie classical Mach-Zehnder 
interferometric approach. Consequently, a length of string proved to 
be a sufficiently accurate measuring device to maintain the two beam 
pathlengths within the coherence length of the laser, a requirement 
for interferogram production. To compensate for the fact that the 
scene beam traversed approximately five inches of glass tunnel walls 


and lucite grids which the reference beam did not, the scene beam 
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was adjusted to be some 2.5 inches shorter than the reference beam. 
Reference beam pathlength was maintained at approximately 138 inches 
throughout the experiment. 

Holograms produced using the basic holographic setup shown in 
Figure 3 exhibited clear, well-defined fringe patterns in nearly every 
instance. In deciding on the final arrangement, several techniques 
were tested to improve upon fringe pattern definition. Horizontal, 
vertical and diagonal translations of the diffuser plate in the scene 
beam were considered, varying from 0.001 inches to 0.005 inches. 

A vertical translation of 0.003 inches yielded clear horizontal fringes 
that were quite easily analyzed. The transverse mode selector aperture 
mas varied from 1.0 mm. to 3.0 mm. in increments of 0.5 mm; best 
lighting of the model resulted with use of a 2.5 mm. aperture. The 
temperature of the cooling water circulated through the laser head 
and etalon was varied from 26.0°C. to 28.0°C. in increments of 

0. oC, ; wath 2TOO. providing the best fringe definition. Finally, 

a variety of beam splitters and lenses were tested prior to final 

Brie ction of the best available optics arrangement for the experiment. 
A 2:1 reference to scene beam strength ratio was found to yield very 
good Relopramee 

Two double exposure holograms were taken for each model view- 
ing angle. The first, labeled a double-static exposure, consisted of 
two exposures ina no-flow condition with a 0.003 inch vertical 


translation of the diffuser plate between exposures. The fringe 
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patterns in this hologram provided a measure of the effect of tunnel 
wall glass, grid plexiglass and model lucite on the subsequent double 
exposure. The second, or static-dynamic, exposure consisted ofa 
no-flow exposure, a 0.003 inch diffuser translation, and finally an 
exposure at flow Mach number 0.937. The fringe deviations recorded 
in the region behind the lucite portion of the model by the double- 
static hologram were measured and subtracted from the fringe shifts 
measured in the corresponding static-dynamic hologram. 

Holograms were produced on Agfa-Gaevert 8E75 holographic 
plates, 4 inches by 5 inches in size. As recommended by Collier, 
et al. [11] the development process included: 

1. Five minutes in Kodak D-19 developer 

2. Thirty seconds ina flowing water bath 

3. Five minutes in standard rapid fixer 

4. Thirty seconds ina flowing water bath 

5. One and one-half minutes in Kodak Hypo Clearing agent 
6. Five minutes ina flowing water bath 

7. Five minutes in methanol bath 

8. One minute ina flowing water bath 


9. Drying 


Be PHOTOGRAPHIC TECHNIQUES 
Normal reconstructions of the original scene were made by 


illuminating the holograms with a seven milliwatt continuous wave 
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helium-neon laser beam at a wavelength of 6328 Angstroms. There 

was some slight distortion in the reconstructed scene because of the 
difference in wavelengths of the original scene beam and the reconstruc- 
tion beam; however, the effect was almost totally negated by shrinkage | 
of the holographic plate emulsion during the development process. 

A common technique of image reconstruction was employed, 
utilizing a conjugate reference beam to reilluminate the exposed 
hologram, as shown in Figure 16. The resulting scene was recorded 
on photographic film. Individual points on the photograph are produced 
by a series of non-parallel rays originating from various source points 
on the diffuser plate in the scene beam. Using a reilluminating beam 
of small diameter has the effect of a small aperture at the focal point 
of the imaging lens, filtering out all but a set of nearly parallel rays, 
as shown in Figure 17. The real images produced in this manner have 
a large depth of field, permitting simultaneous projection on the film 
of front and rear grids, the model and the fringe patterns. The imag- 
ing lens was focused as near to the plane of the fringes as possible, 


producing photographs at various planes of constant Zoe 


Cc. DATA REDUCTION 

Photographic interferograms were obtained using the arrangement 
shown in Figure 16, with the camera viewing screen in the position of 
the realimage. The line of sight in the plane desired was achieved 


by translating and elevating the hologram until common points on the 
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front and rear grids were aligned. The Graphic View camera, with 
the wide angle lens aperture set fully open at f7.8, was adjusted to 
yield the best focus on the fringe plane. Exposure times of from 1/5 
to 3/4 seconds were used to produce workable interferogram photo- 
graphs on Polaroid Type 55 P/N film. 

Fringe shift analysis was accomplished on 8 inch by 10 inch 
enlargements of the 4 inch by 5 inch film used to record the images. 
The enlargements were placed face down ona light table, and the 
fringes, model contours and shock wave were traced on the back 
surface at the desired cross-sectional plane. Fringe shift values 
were recorded by measuring the distance between (1) the intersection 
of the hypothetically undeviated fringe and the cross-sectional plane 
and (2) the intersection of the deviated fringe with the same plane. 
Fringe shifts so obtained were corroborated by placing the negative 
in a photographic enlarger and tracing the lines of interest directly 
onto graph paper. 

The known model fuselage diameter of 1.1 inches was compared 
with that measured in each individual a ere to yield magnification 
factors relating projected dimensions to actual dimensions. These 
factors were then used as corrections to the fringe shift measurements. 
A grid reference point located on the cross-sectional plane of interest 
served as the datum for all fringe shift measurements. A base point 
was located at the intersection of the cross-sectional plane of interest 


and the body longitudinal axis. Fringe shift measurements were 
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corrected using this base point as the new datum so that the inversion 
circle was properly centered on the body axis. The radius of the 
inversion circle was selected to be nearly equal to the semi-span of 
the wing. Fringe shift measurements were converted to fringe 
numbers using the average free stream spacing, while fringe locations 
were nondimensionalized using the inversion circle diameter. From 
the data so obtained, the radial variation of fringe number was plotted 
for each viewing angle. Fringe numbers at 201 equidistant points, as 
required for input into Mode 3 of the computer program, HOLOFER, 
were recorded from the resulting curves. Further details concerning 
this inversion computer program are outlined in Appendix B. A typical 
reduction of an interferogram to obtain the radial variation of fringe 


number at a particular cross-sectional plane is detailed in Appendix A. 
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V. EXPERIMENTAL RESULTS AND DISCUSSION 


A pair of double exposure holograms was taken of the model at 
11% degree intervals through a 180 degree field of view. Experimental 
data from the wind tunnel runs are recorded in Table l. Initial result- 
ing density patterns indicated relatively smooth contours across 
adjacent intervals; the interval was therefore doubled to 2212 degrees to 
simplify and speed the analysis. Fringe data were first inserted into 
the inversion computer program along nine lines of sight in the 180 
degree field of view. A numerical comparison of views from 0 degrees 
to 90 degrees and from 90 degrees to 180 degrees verified to within 
0.20 percent the assumption of a single plane of symmetry in the 
experiment. The fringe data input was then reduced to five lines of 
sight ina 90 degree field of view, as shown in Figure 18. The result- 
ing output was an inverted density field along nine radial lines span- 
ning a 180 degree field of view, with a mirror image on the opposite 
side of the plane of symmetry, as shown in Figure 19.. 

The eine ee photographic interferogram for the 0 degree 
view, along with its corresponding double-static interferogram, is 
shown in Figure 20. The diffraction effects caused by the presence 
of the lucite portions of the model are clearly visible in the double- 
static exposure, where the free stream fringe lines are bent and 


displaced toward the model axis. This displacement was measured 
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and subtracted from subsequent measurements made ongtheystatie— 
dynamic exposure, as outlined in SectionIV.A. Photographic inter- 
ferograms of the remainder of the static-dynamic exposures are shown 
in Figures 21 through 24. Clearly visible and reduceable in nearly all 
views were (1) the region of uniform subsonic flow, commonly called 
the free stream condition, (2) the Ee neition from local subsonic to 
local supersonic flow, and (3) the lambda-type shock wave on the model 
wing. These characteristics are shown in schematic representation 

in Figure 25. 

Contour plots of the density function, as expressed in Equation (8) 
of Section III.A., for successive z-planes of analysis are shown in 
Figures 26 and 27. The cross-sectional plane of analysis for the 
plot of Figure 26 was located at 186.75 mm. from the model nose 
along the longitudinal axis. For Figure 27, the plane of interest was 
195.25 mm. from the model nose. It is apparent from both contour 
plots that the model went to a very small angle of attack under the 
loading forces produced during tunnel operation; this is evidenced by 
the compression of the contour lines above the model and the corres- 
ponding expansion of the contours below the model. Measurements 
made from photographic interferograms confirm this angle of attack 
to be, at most, 0.05 degrees. The closed contours above and adjacent 
to the wing surface in both figures may very well be the result of a 


vortex originating at the intersection of the wing leading edge and the 
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fuselage on either side of the model and traveling aft and outward 
over the wing surface. 

A comprehensive quantitative analysis of the shock wave structure 
was not undertaken, with the exception of estimating the strength of 
the shock by comparison of fringe line separation immediately ahead of 
and aft of the shock wave. Fringe line separation measurements on 
either side of the shock wave were converted first to density informa- 
tion and thence to pressure information, disregarding compressibility 
effects. An approximate strength value of 0.207 was computed using 
the accepted definition of (Pp, - PH) IP 4: This corresponds to a local 
Mach number of 1.08 in the supersonic region just ahead of the shock 
wave. Qualitative shock wave analysis resulted in the construction of 
a three-dimensional structural representation as shown in Figure 28, 
using input information from several interferogram viewing angles. 
While location of the leading and trailing edges of the lambda-type 
shock wave was very accurate, interior structure was largely indefin- 
able due to "smearing"! and blurring of the fringes transiting the shock 
wave 1esiel i. 

As a preliminary step to possible future studies in this field, 
photographic interferograms were made from holograms produced 
with the aerodynamic model set at small angles of attack. Orienta- 
tions included angles of attack of five and ten degrees, with roll 
angles varying from zero to ninety degrees. Although the holograms 


themselves were of very good quality, the photographic reproductions 
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were relatively poor due to the fa¢t that an inferior photographic 
arrangement had to be used. They were therefore omitted from this 
report. It was inferred from the holograms, however, that a complete 
study at angle of attack using the basic procedures followed in the 
present study would be both totally feasible and rewardingly fruitful. 
The original character of the experimental data prevented com- 
parison with published results. Qualitative studies of transonic 
phenomena are widely available, and the general characteristics of 
the resulting density field and shock wave structure serve to bear 
out the schematics based on theoretical and mathematical models. 
Moreover, the self-testing mode of the inversion computer program, 
HOLOFER, verified the consistency and reproductibility of the result- 
ing density distributions to within 2.0 percent through proper choice 
of the input parameters, primarily the slope-matching parameter &. 
The errors encountered in the final results are due primarily to 
errors in the fringe data input to the inversion program. The 
intrinsic presence of laser speckle, the extended pathlengths of the 
Beene and reference beams and the unavoidable beam scattering and 
diffraction within the lucite model sections, created difficulty in 
obtaining precisely the slope of the fringe lines behind the lucite 
sections. Fringe spacing measurements in the free stream flow were 
conservatively judged accurate to within 0.5 mm. This assumption 
was quite reasonable since all measurements were effected witha 


scale graduated at half millimeter intervals. This figure of 0.5 mm., 
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combined with the mean free stream fringe spacing of 5.197 mm. for 
all interferograms, indicated measurement accuracy to within one 
tenth (0.1) of a fringe. The mean systematic error of the free stream 
spacing in each view was computed to be a maximum of 3.9 percent. 
Associated with this systematic error was a random error of 2.1 per- 
cent in the measurement of fringe shifts in each view to a conservative 
accuracy of 0.5 mm. The resulting error for each viewing angle was 
therefore a maximum of 6.0 percent, found by merely adding the two 
types of error for each view. The minimum error limit was found by 
considering the error resulting from the reproduction of the same 
interferogram view five separate times. Statistically, with 6.0 per- 
cent error in each view, the composite error for the repeated view is 
2.6 percent. As five different views, or lines of sight, were used for 
the data between zero and ninety degrees, the final total error in the 
analysis was therefore in the interval between 2.6 percent and 6.0 
percent. To insure contour clarity and guard against overlapping, 
the maximum error figure of 6.0 percent was used in construction of 
the plots shown in Figures 26 and 27. In general, the rather large 
fringe shifts led toa ae low mean fringe sensitivity value of 0.1259. 
This coefficient indicated a resulting density function (Equation (8)) 
inaccuracy of less than 1.5 percent for a fringe shift measurement 
misreading of 0.5 mm. 

Physical limitations of beam diameter and hologram plate area 


dictated the choice of an inversion circle diameter somewhat smaller 


28 


ear! pu as ret 2 Yo" seionde “otal 4 
Pe a ; 
re : with is oe eaatey 


go 350 sant ab: ro 209TH suiedos ga 
hego'tad €: ton ajonieem « 56 Raga a = “Ms 


on 
ie 


cog 2 5 lo, & ial tent odaw toes vacate 
rier i98LUNd HO wi ian as aNidea ‘oyueta Yer a 
wl its, esliaee dose col vols qulltnaee tei ee 
ie die Vi Pine ed said, Jheathg G bko en 
AG “ew GST POTTS iputiceliadie oat seo ita ea 
nica’ ath Iq nokechdages sits ita ah yaltlnges iatie's 
too 0.8 lw -ntlastelgere <a nit dca ahaa 


wily Qofrsties sri aor out! wiisoqmmas Oily waive 


=o 


pl Sou Ginw dilate Ge aanlt to were seeididetoalabae™ i wha 
% Senge MEF 4a 


aA i hy 


sii nl seer! Tepsy tyunhsads ,o soreeh chink fra aesny 
» ideatein 4 an avec aal ce 
yOorlioe> se6p8 hear Reel na wsttiaainan eunnbs 
* 
Ste it gt ee “7 Lo to oxaglie wera iG 
ol ‘radipx ade. See Ae “6 38 oan Oe 
whe 4 wT wets hla eal ae wor bel ass 
} op ort) ies wh i" eee Lars Range. : 
‘ath elon ee: f wilds 496i Sol eh inekaiadl 
| | ‘ve ham sion cud 
Se) tian AEH! aed eo 0 sei Eeaanyt . 
ne, hs al . a i gevndn siod en aye! a9 he meat i tenants 


3 , ‘ | : ma Ma f i « rh ie ey - ‘ i ; r 
me agian. A ela, OP 


than the full data circle normally used in the finite fringe procedure. 
This, in effect, introduced an inconsistency in reference density into 
the analysis since the density on the selected circle and immediately 
outside it was not the calculated Qo ; the density function was there- 
fore not zero outside the actual inversion region. To alleviate this 
inconsistency, a new, updated persone density was computed for 
each cross-sectional plane of analysis by averaging the density values 
on the selected circle from the first inversion process. The actual 
reference densities used were P = 12 %71 me/cce ter the 186.75 mm. 
plane and ve —HlG42 ump) ce torithedl95.25 mam. plane, his pro- 
cedure was justified since all density values between the selected 
circle and the full data circle were constant to within approximately 
fifteen percent. The updated reference densities were then used to 
produce the final output density field. The net effect was a scaled, 
uniform shift toward density function values slightly lower than those 


computed on the basis of the original reference density. 
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VI. CONCLUSIONS 


The finite fringe procedure for the production of holographic 
interferograms has been applied successfully to the determination of 
the three-dimensional density distribution of the flow near the wing- 
fuselage junction of a partially transparent aerodynamic model in the 
transonic regime. Density contours accurate to within six percent 
enabled a thorough analysis of the flow field to be conducted, high- 
lighting flow characteristics and the presence of the shock wave. 
Subsequent studies of similar models at angle of attack have been 
shown to be entirely feasible. Procedures used in the experiment 
also exhibit promise for the direct analysis of duct and inlet flows as 
well as comprehensive study of shock wave structure. 

The inversion computer program, HOLOFER, was found to be 
adequate in handling a general asymmetric flow field analysis. How- 
ever, it was considered quite cumbersome and difficult to modify for 
various experimental situations. Subsequent analysts will find the 
procedures advocated by Sweeney and Vest [12] for the recording and 
analyzing of interferograms of considerable interest. In addition, the 
efforts of Van Houton [13], who utilized the method proposed by 
Junginger and van Haeringen [14], may prove valuable in reducing 


computer time significantly. 
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FIGURE 2. BECKMAN 210 DATA RECORDING SYSTEM 
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APPENDIX A 


REDUCTION OF AN INTERFEROGRAM 
TO OBTAIN FRINGE SHIFT DATA 

The reduction of data for a cross-sectional plane of interest 
involved a complete analysis of both doubly exposed holograms, and 
their corresponding interferograms, for each viewing angle. 

For each view, or line of sight, the double-static exposure was 
first placed ae down ona light table. The average fringe spacing 
was recorded on the back by measuring the perpendicular distance 
between two sufficiently separated fringe lines in the free stream and 
dividing by the number of spacings in between. The method is portray- 
ed in Figure 29. Of primary importance was the measurement of the 
uniform shift of fringe lines due to the beam diffractive quality of the 
lucite sections of the model; this was taken to be the average distance 
between hypothetically unaltered fringe lines and the corresponding 
shifted fringes at their intersection with the z-plane. 

The static-dynamic exposure was then placed face down on the 
light table and the fringe and model contours traced on the back, as 
shown in Figure 30. Again, the average free stream fringe line 
spacing was measured and checked against the value from the double- 
static exposure. A reference point, or datum, was selected at the 
intersection of a well-defined horizontal grid line and the z-plane 


of interest. Fringe shifts were computed in the following manner: 
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1. The distance from datum to a hypothetically undeviated fringe 
line at its intersection with the z-plane was measured in 
millimeters. 

2. The distance from datum to the actual deviated, or shifted, 
fringe at its intersection with the z-plane was measured in 
millimeters. Measured shifts should be made perpendicular 
to the fringe direction. The present procedure is convenient 
and only introduces an error of about 1% in the overall level 
of the density field. 

3. The raw fringe shift distance was then the distance in 1. above 
minus the distance in 2. 

4. To correct for diffraction error, the uniform shift measured in 
the double-static exposure was then subtracted from the 
distance in 3. 

Fringe numbers were calculated by dividing the shift for each fringe 
by the average free stream spacing for the exposure. Magnification 
factors were computed for each exposure by comparing photographically 
recorded model diameter with actual model diameter. 

Since the datum location varied slightly from exposure to exposure, 
it was necessary to reference all measurements to a central point for 
each plane of analysis. This point was taken to be the intersection of the 
longitudinal axis of the model fuselage and the cross-sectional plane of 
interest. The intersections of shifted fringe lines with the z-plane were 
referenced to this fixed point and demagnified. The resulting distances 
were then nondimensionalized using the selected inversion circle radius. 
Table II contains typical data recorded for the zero degree view at the 


186.75 mm. plane of interest. 
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A plot of fringe numbers versus nondimensionalized fringe | 
location was produced for each viewing angle using data obtained from 
the interferograms. Fringe numbers at 201 equidistant points, as 
required for input tae Mode 3 of the inversion computer program, 
were recorded from the resulting curves. The curve plotted using 


the data from Table II is shown in Figure 31. 
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APPENDIX.B 
APPLICATION OF COMPUTER PROGRAM "HOLOFER" 


The computer program used in this study is designed to invert an 
array of fringe numbers across a field to obtain the associated density 
field using a special adaptation of the inversion technique first proposed 
by C. D. Maldonado [5,6]. Three different modes of operation are 
available to the operator, as described below: 

(a) Mode 1 

This mode provides the basic self-testing capability of the 
computer program. It can either generate its own input density field 
using Subroutine FUNCT or read ina density field through Subroutine 
FREAD. The fringe number array corresponding to the input density 
data is first generated; this information is then used as the input for 
the inversion to obtain the original density distribution once again. 

This self-testing procedure was utilized in the present investigation 
to Aeverming the best value of the scale factor, C{ , required to assure 
accurate density in the region of inversion. 

(b) Mode 2 

This mode generates the fringe array at regular intervals 
across the test field from irregularly spaced fringe values read in 


through Subroutine SHEET. Simulated fringe arrays may be generated 
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iby one of the functions specified in Subroutine FUNCT if NCODE = 1 
is specified. Mode 2 has not been utilized in the present investigation. 
(c) Mode 3 

Mode 3 operation reads ina regularly spaced array of input 
density values directly utilizing Subroutine READ, which is first called 
by Subroutine GARRAY. The parameters in the first two cards preced- 
ing the input fringe data serve to identify the size, location and sym- 
metry of the fringe field. The following parameters were used in 


calculating the asymmetric density field in the present experiment: 


PARAMETER INPUL 
NOF Run Number 
IMAX 201 
JMAX 20 
ISYM i) 

JSYM 1 
IMS | 201 
JMS 5 
Z 0:0, 1-0 
XO 0.0 
Y@ Oro 
PHIS YM 0.0 


Amplifying details and applications of the inversion computer 
program are found in References [3,7,9]. A listing of the program is 


included in this appendix for reference. 
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